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Abstract 

A new fit of a chirally motivated coupled-channels model for meson-baryon 
interactions is presented including the recent SIDDHARTA data on the Is 
^ I level characteristics of kaonic hydrogen. The kaon-nucleon amplitudes gen- 

pH I erated by the model are fully consistent with our earlier studies. We argue 

that a sharp increase of the real part of the in-medium K~p amplitude at 
CNJ ■ subthreshold energies provides a link between the shallow ^-nuclear opti- 

cal potentials obtained microscopically from threshold KN interactions and 
the phenomenological deep ones deduced from kaonic atoms data. The im- 
^ I pact on the A-dependence of the A-hypernuclear formation rates measured 

in reactions with stopped kaons is discussed too. 
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■ 1. Introduction 



The interaction of kaons with nuclear medium is standardly described 
in terms of kaon-nuclear optical potential constructed as a coherent sum of 
kaon interactions with individual nucleons, Vk{p) ~ Fkn P, where Fkn is the 
effective elementary kaon-nucleon scattering amplitude and p stands for the 
nuclear density. For many years the strength of the attractive i^'~-nuclear op- 
tical potential remained a puzzle with two conflicting scenarios. The density 
dependent formulation of the Fk-n amplitude implemented in phenomeno- 
logical fits to kaonic atoms data gave a strong evidence of a deep iir~-nuclear 
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potential [H, in a range Re Vk-{po) ~ —(150-200) MeV at nuclear density 
Po = 0.17 fm~^. On the other hand, models based on a threshold value of the 
KN amplitude calculated within coupled-channels ttTj-KN framework with 
chiral dynamics [2|, jsj, jsf provide an optical potential which is only 
(40-50) MeV deep when nuclear medium effects (Pauli blocking and kaon 
self-energy) are taken into account [6|, In a recent work jsj its authors 
pointed to a strong energy dependence of the chirally motivated in-medium 
K~p amplitude and suggested that the puzzling discrepancy could be re- 
solved by evaluating the amplitude at subthreshold energies rather than at 
the KN threshold. 

In the present paper we demonstrate that the subthreshold energy de- 
pendence of the K~p interaction depends only moderately on particularities 
of the underlying chiral model. Specifically, we show that our new fits of 
model parameters to the fresh experimental data on kaonic hydrogen char- 
acteristics generate KN amplitudes that are quite consistent with those 



obtained with a model [10| fitted to the older kaonic hydrogen data [11|. We 
also present the vrS mass spectrum generated by the model that appears 
consistent with observations of the A(1405) resonance peaking at energies 
around 1400 MeV, rather then around 1420 as advocated by some other 



authors [12|, [13j. Finally, we briefly review the impact of energy depen- 
dence implemented in evaluation of K~N branching ratios on the calculated 
A-hypernuclear formation rates in the {K~^^^^^tt~) reactions IJ], 15 . 



2. Separable potentials model 

The synergy of chiral perturbation theory and coupled channel T-matrix 
re-summation techniques provides successful description of KN interactions 
at low energies. In our approach we employ chirally motivated coupled- 
channel s-wave potentials that are taken in a separable form. 



with a/s, p and p' denoting the total meson-baryon center-of-mass energy, 
initial and final state meson momenta, respectively. Further, Ei, Mi and 
Ui stand for the baryon energy, baryon mass and meson energy in the cm. 
system of channel i. The coupling matrix Cjj is determined by chiral SU(3) 
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symmetry. The parameters fi represent the pseudo-scalar meson decay con- 
stants and the Yamaguchi type form factors gi{p) = + are de- 
termined by inverse range parameters a^. The indices i and j run over the 
meson-baryon coupled channels ttA, ttS, KN, rjA, VjH and i^S, including all 



their appropriate charge states. Details of the model are given in Ref. [10 
where a common value of /j = / ~ 100 MeV was adopted in all channels 
and used as a free parameter in fits to experimental data. Here we follow the 



example of Refs. [16|| and [17| and allow for three different couplings /^r, fx 
and fn- 

The chiral symmetry of meson-baryon interactions is reflected in the 
structure of the Cij coefficients derived directly from the Lagrangian. In prac- 
tice, one often considers only the leading order Tomozawa- Weinberg (TW) 
interaction with energy dependence given by 

a,(v^) = -Cj;^{2^s - M, - M,)/4 , (2) 

with Cj^ standing for the SU(3) Clebsh-Gordan coefficients. The exact 
content of the matrix elements up to second order in the meson cm. kinetic 



energies was already specified in Ref. |2| and followed in Ref. [10|. In the 
present work we adopt an alternate formulation of the next-to-leading order 
(NLO) terms for the s-type (direct) and u-type (crossed) Born amplitudes. 



Following the prescription for these two terms by N. Fettes [18[, [19[ and 
keeping only the contributions up to second order in meson momenta the 
two NLO terms read 

Cffi^ = 

Cfi-'(v/5) ^ CS"^(-p?-p|4£g-). (3) 

where the coefficients C[^'^ are the same as those specified in the Table 11 of 
Ref. [lo|, Mo stands for the baryon mass in the chiral limit and rrii denote 
the meson masses. The main advantages of the adopted formulation are a 
scalar character of the baryon propagator and its angular independence in 
case of the u-term. The vanishing of the s-term appears due to averaging 
over angles in the s-wave and would not occur if our model involved higher 
angular momenta. The resulting form of our u-term requires additional com- 
ments. In a proper NLO expansion, there should be a product of meson 
energies Ui and Uj in the denominator of the last contribution instead of the 
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masses mj and rrij. However, this would make the u-term divergent for very 
deep subthreshold energies where pf ~ —^f- This would lead to unphys- 
ical divergences in our meson-baryon interaction potential ([1]) that should 
be cured by higher than NLO terms. Since we terminate the chiral expan- 
sion at the second order we have to regularize the u-term in some other way 
which we do by approximating the meson energies there by their masses. We 
have checked that this approximation is perfectly justified in a large energy 
interval around the KN threshold, without any observable impact on the 
resultinlg KN amplitudes as far as to vrE threshold when extrapolating to 
KN subthreshold energies. 

When the separable potentials ([T]) are used in coupled channel Lippman- 
Schwinger equation the resulting scattering amplitudes are also of a separable 
form given explicitly by 

(4) 

Here the meson-baryon propagator G{^/s) is diagonal in the channel indices i 
and j. When the elementary KN system is submerged in the nuclear medium 
one has to consider Pauli blocking and self-energies (SE) generated by the 
interactions of mesons and baryons with the medium. Thus, the propagator 
G{^/s) and the amplitudes Fij become dependent on the nuclear density p. 
The intermediate state Green's function is calculated as 

where Pi is the on-shell cm. momentum in channel i and the integration 
domain fij(p) is limited by the Pauli principle in the KN channels. Included 
in the denominator of the Green's function is the sum Hj of meson and 
baryon self-energies in channel i. Since the kaon SE is constructed from 
the resulting KN amplitudes a selfconsistent procedure is required as first 
suggested by Lutz 20|. In our calculation, following Ref. 0], the baryon and 



pion self-energies were approximated by momentum independent potentials 
V = Vq p/po with real and imaginary parts of Vq chosen consistently from 
mean-field potentials used in nuclear structure calculations and in scattering 
calculations, respectively. Specifically, we adopted Vq = (30 — ilO) MeV, 
V^A = (_3o _ iio) MeV, V^o^ = (30 - ilO) MeV and Vq^ = (-60 - ilO) MeV. 
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2.1. Fits to experimental data 

The free parameters of the separable-interaction chiral models considered 



in Ref. [10[ and in the present work were fitted to the available experimental 



data on low energy KN interactions consisting of 

• the K~p cross sections for the elastic scattering and reactions (see 



references collected in j2|); following the procedure adopted in [lO| we 
consider only the data points at the kaon laboratory momenta Plab = 
110 MeV (for the K~p, K^n, vr+S", 7r~S"'" final states) and at plab = 
200 MeV (for the same four channels plus vr'^A and 7r°S°) 

the K~p threshold branching ratios, standardly denoted as 7, Rc, and 



Rn. 21 



the kaonic hydrogen characteristics, the strong interaction shift of the 
Is energy level AEu and the decay width of the Is level Fi^ provided 
by the recent SIDDHARTA measurement joj; the older DEAR data 



111 were used in Ref. 10 



In general, the chirally motivated models have no problem with repro- 
duction of the low energy K~p cross sections, mostly due to relatively large 
error bars of the experimental data. This goes in line with inclusion of the 
data taken only at two representative kaon momenta to fix the cross section 
magnitude, which turns out sufficient for a good reproduction of the exper- 



imental cross sections in the whole low energy region [10[. The threshold 
branching ratios are determined with much better precision and provide a 
sterner test for any quantitative usage of the models. Another stringent test 



is provided by the new SIDDHARTA measurement [9| of the kaonic hydro- 



gen data which appears consistent with the other older data [2l| on the K~p 
reactions at threshold and low energies. 

In the current work we present results of two fits to the experimental 
data, one performed with only the leading order (LO) Tomozawa- Weinberg 
interaction ([2]) which we denote TWl, and another one that includes the 
other first order plus second order corrections as well and will be referred 
to as NLO30. Since we fit only 15 experimental data (10 cross sections, 3 
branching ratios, 2 kaonic hydrogen characteristics) it is essential to reduce 
the number of free parameters as much as possible. Our TWl fit was obtained 
by varying only one inverse range parameter common to all channels, = 
axw = 701 MeV, and one meson decay constant, fi = frw = 113 MeV. For 
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the NLO30 fit we fixed the couphngs /j at their physical values /a- = 92.4 



MeV, fx = 110.0 and = 118.8 MeV (see references cited in [ITj), and 
the inverse ranges of the channels closed at the KN threshold were set to 
OryA = ttjyso = ctKs = 700 MeV. The remaining three inverse ranges a^A, 
and a^AT were varied together with four low energy constants (denoted as 
(i's, see Ref. [l3] for their specification) of the second order chiral Lagrangian. 
The remaining low energy constants (couplings of the second order chiral 
Lagrangian, parameters bo = 0.064 GeV~^, bp = —0.209 GeV~^) were fixed 
to satisfy the Gell-Mann formulas for baryon mass splittings, to give the pion- 
nucleon sigma term a.j,N = 30 MeV (parameter bo = —0.190 GeV~^) and to 
reproduce the semileptonic hyperon decays (parameters D = 0.80, F = 0.46) 
10|. Thus, the NLO30 model has only 7 free parameters fitted to the data 



while 10 free parameters were_used in our earlier work [10| and even more of 

we 



them were employed in Ref. [17|. Finally, we note that unlike in Ref. [10 
keep the energy dependence of the TW term in the same form as in our TWl 
fit, Eq. ([2]). This correction as well as a use of empirical PS meson decay 
constants introduce additional chiral symmetry breaking effects that can be 
viewed as renormalization of the pertinent chiral Lagrangian quantities, thus 
reducing impact of higher orders in the chiral expansion. In Ref. ^ our TWl 
model was already discussed and used in analysis of kaonic atoms and kaon- 
nuclear quasi-bound states alongside with an older NLO model CS30 taken 
from Ref. l^l . Since the CS30 model was fitted to the older DEAR data 11 



on kaonic hydrogen it is important to check how much are our conclusions 
sensitive to the mentioned modifications of the NLO model and to the new 
experimental data from SIDDHARTA 

In Table [1] we show how our TWl and NLO30 models compare with 
other models in reproduction of the KN threshold data. The first three lines 
are reserved for models that implement only the TW interaction, the next 
four lines for representative examples of models that incorporate all LO and 



NLO terms. Only the present and the models from Ref. [ITj use the new 
SIDDHARTA measurement of kaonic hydrogen characteristics while the 
CS30 [i3 and BNW 3 models were fitted to the DEAR data [nj. The 
model of Ref. ^ was fitted only to the threshold branching ratios and to the 
shape of the A(1405) resonance in the vrS mass spectrum. It can be seen 
that the threshold data are quite well reproduced with the models that are 
based only on the TW interaction. When we incorporate the NLO terms 
the total P^r data point drops significantly from 3.6 to 0.61 for the TWl 
and NLO30 models, respectively. The CS30 model gives /N = 0.78 for 
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Table 1: K^p threshold data calculated in several LO and LO+NLO coupled-channel 
chiral models. The columns show the kaonic hydrogen Is level shift AE'i^ and width Tis 
(in eV; the values marked by an asterisk were derived from the K^p scattering length by 
means of the modified Deser-Trueman formula 22]), and the K^p threshold branching 
ratios 7, i?c, Rn- The last two columns list the / = S'-matrix pole positions zi,Z2 (in 
MeV) related to the A(1405) resonance. The last two lines show the experimental data 
and their errors. 



model 


AEu 




7 


Rc 


Rn 






Z2 




TWl 


323 


659 


2.36 


0.636 


0.183 


(1371,- 


-54) 


(1433,- 


-25) 


TW2 [3] 


275* 


586* 


2.30 


0.618 


0.257 


(1389,- 


-64) 


(1427,- 


-17) 


TW3 [17] 


373* 


495* 


2.36 


0.66 


0.20 


(1384,- 


-90) 


(1422,- 


-16) 


NLO30 


310 


607 


2.37 


0.660 


0.191 


(1355,- 


-86) 


(1418,- 


-44) 


CS30 [10] 


260 


692 


2.37 


0.655 


0.188 


(1398,- 


-51) 


(1441,- 


-76) 


BNW [^ 


236* 


580* 


2.35 


0.653 


0.194 


(1408,- 


-37) 


(1449,- 


106) 


IHW [17] 


306* 


591* 


2.37 


0.66 


0.19 


(1381,- 


-81) 


(1424,- 


-26) 


exp. 


283 


541 


2.36 


0.664 


0.189 










error (±) 


42 


111 


0.04 


0.011 


0.015 











Table 2: The inverse ranges (in MeV) and low energy constants do, do, dp and di (in 
GeV~^) fitted to the available K~p data at threshold and low energies. 

model OttA «7rS ctRN dr, dp di 

CS30 291 601 639 -0.450 0.026 -0.601 0.235 
NLO30 297 491 700 -0.812 0.288 -0.737 -0.016 



the data set which includes the new SIDDHARTA measurement. When the 
parameters of the CS30 model were fitted to the DEAR data HI] the resulting 
/N = 1.3 was almost twice as large. This difference demonstrates that 
the new kaonic hydrogen data joj from the SIDDHARTA collaboration are 
much more consistent with the other low energy data on K~p scattering and 
reactions. The two parameter sets of our CS30 and NLO30 models (both of 
them giving a.„N = 30 MeV) are compared in the Table [21 As expected, the 
alternate treatment of the NLO terms and use of empirical couplings /, has 
significant impact on the (i-couplings while the inverse ranges are varied 
only moderately. We also note that our d-couplings are somewhat larger than 



those reported for the IHW model [ITJ, though there is no simple relation 



between their and our NLO low energy constants. 
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The separate treatment of the inverse ranges in different channels and a 
relatively small value of a^A in our NLO models may require a comment. 
First of all we note that only a small part of the improvement observed when 
going from the LO to the NLO30 fit is due to the separate treatment of the 
inverse ranges in the NLO30 model. In a complementary NLO fit with fi 
fixed at physical couplings and only one value of inverse range used in all 
channels and varied together with the NLO d-couplings the fit to the data 
gave x^/^ = 1-1 for c^i = c^NLO = 481 MeV. Although the quality of this fit is 
quite good in general it is lacking in respect to the resulting kaonic hydrogen 
Is-level width Lis = 763 eV, which is by about two standard deviations 
larger than the SIDDHARTA value (the energy shift AEis = 275 eV would 
be fine). The independent variation of inverse ranges in channels open at the 
KN threshold allows to remedy this nuisance. It is also a common practice 
in other KN chiral models jif, Q, 17| to use channel dependent subtraction 
constants to regularize the intermediate state integral. The large negative 
values of the vrA subtraction constant reported in the NLO fits of Refs. j3] 
and [13] are equivalent to a very small cutoff momenta 23 1 which in turn 
implies small values of inverse ranges when the intermediate state integral 
is regularized by means of the Yamaguchi formfactors. Thus, although the 
relation of our inverse ranges to the subtraction constants is highly nontrivial, 
the small value of a^A appears in qualitative agreement with the values of 
the pertinent subtraction constant obtained in other NLO fits of the K~p 
data. 

Finally, we would like to emphasize the importance of calculating the 
kaonic hydrogen characteristics directly by solving the kaon-proton bound 
state problem rather then deriving them from the K~p elastic amplitude 
extrapolated to the threshold. Standardly, the kaonic hydrogen shift AEis 
and width Fi^ are derived from the K~p scattering length a^-p by means of 



\ p 

a modified Deser-Trueman formula 221 taken to the second order in 



-p- 



However, we have shown [5[ that the precision of the relation is on a level 
of 10%. Since the experimental precision of kaonic hydrogen characteristics 
is getting close to this level one should either expand the relation reported 
in Ref. [22| to higher orders of a^-p or solve the K~p bound state problem 
numerically as we do here and in Refs. [sj and 10 



Also listed in Tabled] are the positions zi,Z2 of the two / = S*- matrix 
poles that reside on the second Riemann sheet [— , +] of the complex en- 
ergy manifold, where the signs are those of the imaginary parts of the cm. 
momenta in the vrS and KN channels, respectively. It is remarkable that 
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all the TW models listed in the table are in a relatively close agreement on 
the position of the upper pole Z2. This agreement is spoiled, at least in our 
model (unlike in Ref. 13]), when the NLO corrections are included in the 



inter-channel couplings. In contrast, the position of the lower pole zi exhibits 
model dependence already in the TW models. Generally, it is located much 
further away from the real axis than the pole Z2- The pole Z2 is usually rele- 
gated to the subthreshold behavior of the K~p amplitude and to the A (1405) 
resonance observed in the vrE mass spectrum in KN initiated reactions. 
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Figure 1: The ttE mass distribution. Our results are compared with the experimental 
data taken from Refs. [131 j [25] and 26] with the experimental bars at each energy shown 
in this order. See the text for explanation on the theoretical curves obtained with the 
NLO30 model. 

The Figured] visualizes the ttE mass distribution computed for the NLO30 
model. The solid curve was obtained by tuning the relative couplings of the 
vrS and the KN channels to an / = source to get a peak at 1395 MeV, see 



also Refs. [23[ and [10|| for details. Just for a reference we also show the spec- 
tra obtained by assuming that the / = resonance originates only from the 
ttS channels (dashed line in Fig. [T]) or that it is formed exclusively from the 
KN channels (dotted line). These two lines represent a kind of boundaries 
on the shape and peak position of the spectra in a situation when the low 
energy constants are fixed at the values obtained in our NLO30 fit. The ex- 
)erimental data shown in Figure [1] come from three different measurements 
all exhibiting a prominent structure around 1400 MeV. As 



25, 26 



the observed spectra are not normalized we have rescaled the original data 
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as well as our computed distributions to give 1000 events in the chosen en- 
ergy interval (from 1330 to 1440 MeV). The three measurements give ttE 
distributions that look mutually compatible. However, as the experimental 
spectra contain admixtures of / = 1 and 1 = 2 contributions they cannot 
be compared simply with our theoretical predictions based on an idea of the 
/ = source. In reality, a meaningful comparison of the measured ttE spec- 
tra with theory would also require a proper treatment of the dynamics and 
kinematics of the particular reaction, i.e. an application of an appropriate 
reaction model. Since this is out of scope of the present work the comparison 
of the theoretical and experimental data is presented in Figure [1] only for 
illustration. We also mention that the K~p "E^ti^tt^ data measured by the 



Crystal Ball Collaboration [27[, which are not shown in the figure, yield a 
slightly different distribution with a peak structure around 1420 MeV. The 
two identical pions in the final state of this reaction complicate a direct com- 
parison with the other experiments, so we also find it questionable to relate 



our computed line-shape to the one observed in [27 



2.2. Free space and in-medium KN amplitudes 

The energy dependence of the K~N amplitudes in vacuum and in nu- 
clear medium was already discussed extensively in Ref. [sj where the reduced 
amplitudes (stripped off the form factors Qi) were presented. Here we pre- 
fer to show the full on-shell amplitude -F^jy ^^^^1 anticipate purely imaginary 
momenta pf^j^ = i | | for the subthreshold energies. In the present Fig- 
ure [2] we compare the energy dependence of the elastic K~p amplitude in 
the free space as generated by three different models. The pronounced peak 
in Im Fx-p and the change of sign in Re F^-p point to the existence of a 
quasi-bound state related to the A(1405) resonance closely below the K~p 
threshold. This feature is common to all models based on chiral dynam- 
ics which, in combination with couple channel re-summation, generate the 
resonance dynamically. The three models employed in our calculations lead 
to very similar KN amplitudes above the threshold and are in qualitative 
agreement at subthreshold energies as well. Interestingly, the Figure [2] also 
demonstrates that the NLO30 model fitted to the SIDDHARTA data 0] 
leads to K~p amplitudes that are in close agreement with those obtained for 
the CS30 model fitted to the DEAR data [ll|. Especially, the real parts of 
the amplitudes generated by the CS30 and NLO30 models are very close to 
each other when extrapolated as far as to (and even below) the vrS threshold. 
This feature may be explained by recalling that in fact the CS30 model was 
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not able to reproduce the DEAR data and generates kaonic hydrogen char- 
acteristics that are more compatible with the SIDDHARTA results. In other 
words, the energy dependence of the K~p amplitude is to some extent fixed 
by the very precise threshold branching ratios and by the low energy scatter- 
ing and reaction data (see Ref. 28| for a detailed analysis). The significance 
of the new SIDDHARTA data can be judged in terms of putting additional 
constrains on the models and reducing the theoretical uncertainties when 
extrapolating the KN interaction to subthreshold energies 17 . 
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Figure 2: Energy dependence of the real (left panel) and imaginary (right panel) parts 
of the elastic K~p amplitude in the free space. Dashed curves: TWl model, dot-dashed: 
CS30 model, solid curves: NLO30 model. 

Another point worth mentioning is related to the positions of the poles 
of the S-matrix shown in the Table [TJ Standardly, the threshold behavior 
of the K~p amplitude is viewed as strongly affected by the higher of the 
/ = poles, here labeled as Z2- We note that this pole is located at rather 
varied positions for the three different models depicted in Figure [2J For our 
CS30 model the Z2 pole is found even above the KN threshold and further 
from the real axis, so it might be the lower zi pole that affects more the 
subthreshold K~p interaction in this case. Since all our models lead to very 
similar K~p amplitudes at energies at and above the threshold, it looks that 
the positions of the / = poles cannot be unambiguously determined from 
the current experimental data. In this respect our observation is similar 
to the one made by Shevchenko 
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who arrived at a conclusion that the 
available experimental data can be described equally well by models that 
incorporate either one or two / = poles. 

In Figure |3] we show the energy dependence of the elastic K~p amplitude 
in nuclear medium for the nuclear density p = pq. Both the Pauli blocking 
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Figure 3: Energy dependence of the real (left panel) and imaginary (right panel) parts 
of the elastic K^p amplitude in nuclear medium at the nuclear density p = pQ. Dashed 
curves: TWl model, dot-dashed: CS30 model, solid curves: NLO30 model. 



as well as the meson and baryon self-energies were included when computing 
the amplitudes in a selfconsistent way. Standardly, only 5-10 iterations are 
required to achieve the selfconsistency The qualitative behavior of the 
amplitudes is once again independent of the model used in the calculations. 
It is remarkable that the selfconsistent treatment leads to even smaller differ- 
ences between the CS30 and NLO30 models than they were in the free space, 
especially in the subthreshold energy region. However, the prominent sub- 
threshold increase of in-medium KN attraction was not observed in earlier 
in-medium calculations \d\ which get substantially different results than ours 
already when only Pauli blocking is accounted for. At the moment we have 
no real explanation for the difference, maybe the on-shell treatment of the 
intermediate state propagator used in Ref. [6| is not so well justified when the 
system is submerged in nuclear medium. On the other hand, we were 
able to reconstruct fully the results of Ref. [30|] (obtained without accounting 
for the hadron self-energies) when we switched to their parameter set. 

As it was already demonstrated in Ref. the in-medium dynamics of 
the A (1405) resonance is responsible for the rapid increase of the real part 
of our K~p amplitude at energies about 30 MeV below the KN threshold. 
While the Pauli blocking pushes the resonance structure above the threshold 
the kaon self-energy is responsible for moving it back to energies where it 
is located in the free space. However, when the nuclear density is increased 
the relevant / = pole crosses the real axis above the KN threshold, thus 
moving to the [+, — ] Rieman sheet. Since it is now located much further from 
the physical region (on a Rieman sheet not connected with the physical one) 
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the K~p amplitude no longer exhibits the resonance structure characterized 
by real part of the amplitude crossing zero and imaginary part resembling 
the Gaussian shape. 

Finally, we mention that the free-space K~n interaction is weakly at- 
tractive and its in-medium renormalization is rather weak and exhibits little 
density dependence We also note that a proper treatment of the KN 
system submerged in nuclear medium requires an introduction of realistic 
in- medium KN momenta that are used in the Yamaguchi form factors gi{p) 
instead of the on-shell momenta assumed in Figure [31 The details can be 
found in Ref. fi]. 



3. A-hypernucIear production 

As mentioned in Section [T] the energy dependence of the in-medium KN 
amplitudes is crucial for a construction of the kaon-nuclear optical potential. 
In many applications the KN interaction in nuclear medium is probed at 
energies below the KN threshold. The energy shift from threshold to sub- 
threshold energies was estimated in Ref. [sl and an optical potential based 
on our chirally motivated amplitudes was used there to calculate the bound 
state characteristics of kaonic atoms and kaon-nuclear quasi-bound states. 
The energy dependence of the K~p amplitude shown in Fig. [3] leads to an 
optical potential that is shallow at the KN threshold but becomes much 
deeper at subthreshold energies relevant for the i^~-nuclear bound state sys- 
tems. Thus, the shift to subthreshold energies brings the constructed optical 
potential in agreement with phenomenological analysis of kaonic atoms data 
that favors deep optical potential [if. In this section we demonstrate an- 
other effect of the in-medium energy dependence that stems from the energy 
dependence of the K~N branching ratios. 

The new experimental data [1^] on the A-hypernuclear production in 
(ii'g^Qp,7r~) reactions allow to study the A-dependence of the formation rates 
for the p-shell nuclear targets. The calculated capture rates significantly un- 
derestimate the measured ones, the deeper the K~ potential, the smaller is 



the capture rate [31[. However, one can look at relative rates and the A- 
dependence of the production rates where the impact of both the theoretical 
ambiguities and the experimental systematic errors should not obscure our 
observations [l^. Within a framework of the distorted wave impulse ap- 
proximation (DWIA) the nuclear capture rate per stopped kaon Rfi/K can 
be expressed as a product of three terms, a kinematic factor, the elementary 



13 



0.0 



8 10 12 14 
nuclear mass number 



16 



0.0 



8 10 12 14 
nuclear mass number 



16 



Figure 4: The A-dependence of the Isa hypernuclear formation rates, experimental data 
by the FINUDA collaboration . The theoretical rates are normalized to the ""Li experi- 
mental value and were calculated with a phenomenological density dependent kaon-nuclear 
optical potential (dashed lines) and with the chirally motivated ii'^-nuclear optical po- 
tential (solid lines). Left panel: elementary BK{K~N — >■ it A) fixed at the threshold value 
for nuclear density p — po/2, right panel: energy and density dependent BR. 



branching ratio (BR) of the process BR(i^ — > vrA), and a rate per hyperon 
Rfi/Y 3l|- While the BR has standardly been taken as a constant fixed at 



its threshold value, our work demonstrates the importance of considering the 
energy and density dependence of the BR. 

In FigureHlwe show the calculated Isa formation rates (per stopped kaon) 



15| in comparison with the FINUDA data IJ] for nuclear targets from ^Li to 



^^O. The presented theoretical rates were normalized to reproduce exactly the 
experimental rate for the ^Li target, so one can focus on the A-dependence of 
the rates. The calculated rates in the left panel of Fig. H] assume a constant 
BR generated by the CS30 model, taken at the KN threshold and for an in- 
termediate nuclear density p = Pq/2. In this case the A-dependence is solely 
driven by the rate per hyperon that contains the overlap of the initial and the 
final state wave functions. There, the chirally motivated kaon-nuclear optical 
potential (based on the CS30 amplitudes) does a better job than the phe- 
nomenological one (based of the density dependent amplitudes from Ref. [l|). 
The calculated rates presented in the right panel of the figure incorporate 
in addition the energy and density dependence of the BR that appears due 
to an energy shift from threshold to subthreshold KN energies. Once again 
we see that the implementation of an energy shift from the KN threshold 
to subthreshold energies (where the in-medium KN interaction is effective) 
significantly alters the observed picture. Both the phenomenological and the 
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chirally motivated optical potentials are sufficiently deep at the subthreshold 
energies relevant for the evaluation of the energy dependent BR, so they lead 
to similar Isa formation rates. 



4. Summary 

We have demonstrated that the results of SIDDHARTA experiment 
on kaonic hydrogen are in good agreement with other available data on K~p 
threshold branching ratios and on the low energy K~p cross sections. Our 
new NLO30 model is fully compatible with results of our previous analyses 
and improves the general description of the data. 

The several versions of coupled-channels separable potential models con- 
sidered in our work provide KN amplitudes that exhibit very similar energy 
dependence in the free space as well as in the nuclear medium. Specifically, 
the strong subthreshold energy and density dependence of the K~p ampli- 
tudes, that reflects the dominant effect of the A(1405) resonance, does not 
depend much on a particular version of the model. The observed sharp in- 
crease of K~p in-medium attraction below the KN threshold is a robust 
feature common to all considered models. It is prominent already in the 
TWl model that employs only the leading order TW interaction. The NLO 
contributions present in the CS30 and NLO30 models improve significantly 
the quality of the fit to experimental data but their impact on the resulting 
KN amplitudes is only moderate. 

Since the i^~-nuclear interaction probes subthreshold KN energies where 
the K~p in-medium amplitude exhibits much stronger attraction the result- 
ing i^'~-nuclear optical potential becomes much deeper than when it were 
constructed from the amplitudes taken at the KN threshold. The mecha- 
nism of constructing the optical potential from subthreshold KN energies 
allows to link the shallow ^-nuclear potentials based on the chiral KN am- 
plitude evaluated at threshold and the deep phenomenological optical poten- 
tials obtained in fits to kaonic atoms data. The relevance of this finding to 
an analysis of kaonic atoms and quasi-bound ^-nuclear states was already 
investigated in Ref. |8|]. The results are also discussed in separate reports of 



this journal [32|, [33 



The subthreshold energy and density dependence of the in-medium ele- 
mentary branching ratio BIi{K~N — ?■ ttA) has a significant impact on the 
A-hypernuclear production rates observed in reactions with stopped kaons. 
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The implementation of the effect in our DWIA calculations leads to A de- 
pendence of the formation rates that is in very good agreement with the one 
observed in the FINUDA experiment performed on p-shell nuclear targets. 
It also brings in agreement the results obtained with i^r~-nuclear optical po- 
tentials based on either the chirally motivated or the phenomenological den- 
sity dependent amplitudes. Unfortunately, the magnitude of the computed 
DWIA rates remain much lower than the rates established experimentally. 
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